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Abstract 

We have calculated inclusive one- and two-jet production in photon-photon collisions 
in next-to-leading order superimposing direct, single resolved and double resolved cross 
sections. The results are compared with recent experimental data from the TOPAZ and 
AMY collaborations at TRISTAN. Good agreement is found between experiment and the 
^ , theoretical results. 
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1 Introduction 



Phenomena involving high transverse momenta (pr) in 7-7 collisions with the incoming pho- 
tons being real or quasi-real have been studied for many years now, both experimentally and 
theoretically Emphasis has been on 7^7 collisions at relatively high center of mass energies in 
order to obtain information beyond the low pr region which is determined by soft physics. This 
has been achieved recently at the TRISTAN |l], |2| and LEP j3| storage rings. Two TRISTAN 
collaborations, TOPAZ |TJ and AMY have presented data for inclusive one-jet and two-jet 
cross sections with pt's up to 8 GeV. This is the domain of perturbative QCD calculations 
which are expected to be valid when there is a large energy scale in the interaction. Then the 
lowest order diagram is that of the quark parton model (QPM) (see Fig. la) also referred to as 
the direct contribution in leading order (LO) QCD. Here the two photons couple directly to the 
charge of the bare quark with no further interactions involved, so that the final state consists of 
just two jets. In next-to-leading order (NLO), i.e. 0(a 2 a s ), there exists a contribution as, for 
example, the one in Fig. lb where a gluon is emitted from the internal quark line. This contri- 
bution becomes singular when the momentum of the outgoing quark (antiquark) is collinear to 
the momentum of the upper (lower) photon. This divergent contribution is dealt with by intro- 
ducing the scale dependent photon structure function at the upper or lower vertex. This leads 
to the so-called 'single resolved' (SR) contribution depicted in Fig. lc in LO, being 0(aa s ), 
when the photon structure function is factored out. NLO corrections to the single resolved 
cross section, for example, produce contributions with a gluon emitted from the internal quark 
line in Fig. Id which becomes singular in the collinear limit. These singular pieces are absorbed 
as before in the scale dependent photon structure function now at the lower (upper) vertex. 
This is the 'double resolved' contribution (DR) with photon structure functions at the upper 
and lower vertex, respectively, as shown in Fig. le. In addition to the two high px jets the 
single resolved (double resolved) contributions have one (two) low pt "beam jets" or photon 
remnant jets. This was first discussed in the paper by Brodsky et al. [|J who showed that 
resolved photon cross sections are of the same order as from LO direct contributions, i.e. the 
QPM process. This work has recently be extended by Drees and Godbole ||. 

In all these processes, direct (D), single resolved (SR) and double resolved (DR) there is a 
hard subprocess in which two partons scatter at high px which can be calculated in perturbative 
QCD. The subprocess cross sections are convoluted with luminosity functions that give the 
distributions of the virtual photons in an electron and one or two photon structure functions 
that give the distribution of quarks and gluons inside a photon. The topological structure of 
these cross sections is equivalent to cross sections for jet production in ^p collisions. There the 
direct (resolved) process corresponds to the SR (DR) component in 77 collisions. The proton 
structure function must be replaced by the photon structure function. 

In LO calculations the three cross sections, D, SR and DR are superimposed, and rates 
for topologies with two high p? jets plus one or two photon remnant jets can be predicted. It 
is well known that LO calculations give only qualitative results. In this order the predictions 
depend highly on the factorization scale, in particular for the DR contribution. Furthermore 
LO results do not depend on any jet definition which is present in NLO and in the experimental 
data. Therefore all three processes must be treated in NLO, i.e. with higher order corrections 
included both in the photon structure function and the hard subprocesses in a consistent way. 
Such calculations were done by Aurenche et al. for the inclusive single-jet cross section at 
y/S = 58 GeV f|. They compared their results with the TOPAZ data and found good 
argeement up to the highest pt by slightly adjusting the non-perturbative input of the photon 
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structure function still consistent with the data on K7. Besides the inclusive one-jet data from 
TOPAZ [|J and AMY || there exist also measurements of the inclusive two-jet cross sections 
by the TOPAZ [|TJ and the AMY collaborations. These data have been compared so far 
only to LO calculations by the authors assuming different sets of photon structure functions. 
In this work we make an effort to predict these two-jet cross sections in NLO assuming the 
same kinematical constraints as in the two TRISTAN experiments. Unfortunately we are not 
in the position yet to calculate all three components, the D, SR and DR cross sections up to 
NLO. We have only the NLO formalism for the direct and the single resolved cross section at 
our disposal. Since they yield, depending on p T , the major part of the two-jet cross section at 
58 GeV we shall estimate the DR cross section with a lowest order calculation including some 
k factor. This k factor is estimated with the inclusive one-jet cross section by comparing with 
the NLO prediction. 

In section 2 we explain the formalism used to calculate the inclusive two-jet cross section. 
Section 3 contains the comparison with the TOPAZ and AMY data and a discussion of the 
results. 



2 Inclusive Jet Cross Sections 

As explained in the introduction the inclusive jet cross section at large px consists of three parts 
according to how the incoming photons take part in the hard subprocesses. If both photons 
couple to the quarks in the hard scattering one defines the direct cross section a D (see Fig. la). 
For example, the inclusive two-jet cross section has the following form: 

d 3 v D = ( d 3 a D (^ -> jet! +jet 2 ) \ | K D/ R m \ m\ 

dpTdr\\dr\2 \ dpTdr\\dr\2 /LO ^ 

In ([I]) pt denotes the transverse momentum of the measured or trigger jet with rapidity rji, 7/2 
is the rapidity of another jet such that in the three-jet sample these two jets have the highest 
and second highest px- The first term on the right-hand side of (|l|) stands for the LO cross 
section and the second term is the NLO correction which depends on the factorization scale 
M, at which the initial state collinear singularity is absorbed in the photon structure function. 
The variable R is the usual parameter defining the size of the jets with transverse momentum 
Pt, rapidity 77 and azimuthal angle 0. When two partons fulfill the snowmass constraint |7|] 
with the cone size parameter R they are recombined in one jet. The same jet definition is used 
in the TOPAZ and AMY analysis. 

The NLO corrections K D (R, M) are calculated with the phase-space slicing method using 
for the seperation of the 77 — > 2 jets and 77 — > 3 jets cross sections an invariant mass cut-off 
y, defined as 2piPj < ys, where s is the partonic center of mass energy squared. In the direct 
process the cross section for 77 — > qqg has soft, initial and final state collinear singularities. The 
77 — > QQ9 cross section is integrated over these singular regions up to the cut-off value y which 
isolates the respective singularities which are cancelled against the singular contributions of the 
virtual corrections to 77 — > qq and the subtraction term at the scale M which are absorbed 
into the photon structure function. Outside the cut-off region controlled by the parameter y 
we have genuine qqg final states. In this contribution two of the partons are recombined if 
they obey the snowmass constraint. The LO contribution and the NLO correction inside the y 
cut-off contribute to the two-jet cross section together with the contributions inside the cone 
with radius R. That part of the qqg contribution not fulfilling the cone recombination condition 
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stands for the 3-jet cross section of which we have calculated the inclusive two-jet cross section 
as a function of pt, f]i and r) 2 where p? is the transverse momentum of the measured jet or 
trigger jet with rapidity 971 and r\ 2 is the rapidity of another jet so that both p Tl and px 2 are 
the jets with highest pr- For the two-jet cross section p^ x = pt 2 . This inclusive cross section 
is independent of the invariant mass parameter y used to seperate those regions of phase space 
that contain the soft and collinear singularities. However the inclusive two-jet cross section 
depends on the cone size R. 

If one of the two incoming photons interacts via the quark or the gluon component in its 
structure function we define 

d 3 a SR 



dpT<irjidrj2 



dp T dr]i7] 2 



This cross section has the same structure as the direct inclusive two-jet cross section for pho- 
toproduction ■yp — > jet x + jet 2 + X, where the photon structure function in @ is replaced by 
the proton structure function. This cross section has been calculated recently with the phase 
space slicing method with invariant mass cut slicing ||. The developed formulas for the two- 
jet photoproduction cross section can be transformed easily to the corresponding 7-7 process. 
The first term in each bracket of (0) is the LO cross section and K-L R , K^ R stand for the NLO 
corrections which depend on R and the two scales \x (renormalization) and M (factorization). 
The cross section in (0) has two separate terms, depending whether the first or the second 
photon interacts through its structure function F i j 1 . 
When both photons are resolved we have 

d 3 a DR r r 

Jd Xl Jdx 2 F ih (x 1 ,M)F jh (x 2 ,M) 



dpTdr]idr] 2 




f, -(/') V DH 



2tt 



K^(R,M,») (3) 



This cross section has the same structure as the NLO resolved inclusive two-jet cross section 
for photoproduction yp — > jet x + jet 2 + X, where the photon structure function at one vertex 
is replaced by the proton structure function. For this cross section the higher order corrections 
K® R (R, M, n) have not been calculated yet. However they are available for the inclusive single- 
jet cross section M and are transformed to the 7-7 case. Concerning the inclusive two-jet 
cross section in the double resolved case our strategy is as follows. First we have calculated the 
inclusive one-jet cross section in NLO, i.e. with the terms K® R , and have compared it with 
the LO cross section, so that we know the k factor for the one-jet cross section. The same k 
factor is applied to the two-jet cross section which is calculated in LO only. Since for p T > 4 
GeV the DR contribution is less than 20% this estimate of the two-jet double resolved cross 
section is sufficient. 

It is obvious that the NLO corrections K D , Kf R and Kf R in ([I]), and (|3|) depend on 
the same kinematic variables p?, 771 and r] 2 as the LO cross sections. We specified only the 
dependence on the cone size R, on the factorization scale M and on the renormalization scale 
ft. 



3 



Before we compare our results with the TOPAZ and AMY data we have made some tests 
of the NLO corrections to the one- and two-jet cross sections. First we checked that these 
cross sections are independent of the cut-off y if y is chosen small enough. This was the case 
for y < 10~ 3 in all considered cases. For y above 10~ 3 we observed some small y dependence 
which is caused by our approximation that we neglected contributions 0(y) in the analytical 
contribution to the two-jet cross section. 

Second we tested that, first, the sum of NLO direct and the LO single resolved cross section 
and, second, the sum of the NLO single resolved cross section and the LO double resolved cross 
section are independent of the factorization scale M. This test was performed for the one- and 
two-jet cross section seperately under kinematic conditions as for the TOPAZ data. Similar 



checks were done earlier for the photoproduction one-jet cross section 1 10 



3 Comparison with TOPAZ and AMY Data 

In e + e~ collisions at TRISTAN the jets are produced via the exchange of two quasi-real photons. 
The spectrum of these small Q 2 photons is described by the Weizsacker- Williams approximation 



F j/e (z,E) = - ^^-^ ln r^ 1 -^ (4) 




where O max is the maximally allowed angle of the electron (positron) which experimentally is 
realized by the anti-tagging counters. E is the beam energy of the electron (positron) and 
z = Ery/E is the fraction of electron (positron) energy transferred to the respective photon. 

In this approximation the inclusive two-jet cross section e + + e~ — * e + + e~ +jet 1 + jet 2 + X 
is obtained from 

d 3 a(e + + e~ -> e + + e" + jet x + jet 2 + X) _ 
dpTdrjidr]2 

P xp , rf 3 a( 7 + 7 -> jet! + jet 2 + X) 
d Zl jdz 2 F j/e (z u E)F, /e (z 2 ,E) (5) 

where the cross section on the right-hand side stands for the interaction between two real 
photons. 

Depending on the experimental conditions we encounter specific values for 6 max and even- 
tually for the integration range of z 1 and z 2 , respectively, in (|5|). 

Before presenting the results we must specify the parton distributions in the photon, Fj/ 7 , 
which appear in (|2|) and (|3|). We have chosen the NLO set of Gluck, Reya and Vogt (GRV) 
in the MS scheme |ITJ. In [^TJ the photon structure function is given in the so-called DIS 
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scheme. The transformation to the MS scheme is known and given in ||11|| . This means that 
the direct, single resolved and double resolved cross section must also be calculated with the 
MS subtraction. We have tested that the DIS 1 scheme for the photon structure function of 
GRV leads to the same results. We choose all scales /i = M = pr and calculate a s from the 
two-loop formula with Nf = 4 massless flavours with A^ = 0.2 GeV equal to the A value of 
the NLO GRV photon structure function. The charm quark is treated also as a light flavour 
with the boundary condition that the charm content of the photon structure function vanishes 
for M 2 < m 2 (m c = 1.5 GeV). 

The TOPAZ data are obtained for O max = 3.2° and z\, z 2 < 0.75 PJ. The rapidities of the 
two jets are restricted to \rji\, \i] 2 \ < 0.7 and R — 1. The same input is used for the calculations 
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of the NLO cross sections. Our result for the inclusive single jet cross sections da '(one-jet) / dpx 
is shown as a function of px in Fig. 2a. We show four curves, the direct, single resolved, double 
resolved and the sum of all three components. All three cross sections are calculated up to 
NLO and the curve for the sum should be compared with the TOPAZ points. The agreement 
is good at small px- At larger px, where the theory is more trustworthy, the theoretical curve 
falls off stronger with increasing px than the data. The error bars on the data include both 
statistical and systematic uncertainties. The effects of the NLO corrections are as follows. If 
we compare with the LO result, where only the NLO corrections in the hard scattering are 
removed, i.e. the LO cross section is calculated with the same NLO GRV photon structure 
function as in the NLO computation we encounter the following k factors: k ~ 0.9 (direct), 
k ~ 1.1 (single resolved), k ~ 1.9 (double resolved) and k ~ 1.0 (sum)[|. We emphasize that the 
k factor depends on the way the LO cross section is evaluated. If we use in LO the one-loop a s 
with the same A value the k factors are: k ~ 0.9 (D), k ~ 0.9 (SR), k ~ 1.2 (DR) and k ~ 0.9 
(sum). In Fig. 2b we have plotted the two-jet cross section G?cr(two-jet)/<ipr as a function of 
Pt- We show again the direct, single resolved, double resolved cross section and the sum of 
all three contributions. Over the whole range of px the agreement with the TOPAZ data is 
excellent. We notice that the double resolved cross section falls off with px somewhat stronger 
in relation to the direct and single resolved components as in the one-jet cross section. The 
double resolved cross section is calculated in LO with the same fc-factor as obtained in the 
inclusive one-jet cross section. For the direct and single resolved cross section the k factors for 
the one- and two-jet cross section are approximately equal, so that we can expect this also for 
the double resolved cross section. 

The same plots but now with the kinematical contraints of the AMY experimental data, 
©max = 13°, Zi, z 2 < 1 and \rj\ < 1 for the one-jet cross section and |r^i|, |r/ 2 | < 1 for the 
two-jet cross section are shown in Fig. 3a,b. All other input is as for the TOPAZ curves. The 
k factors are the same as for the TOPAZ cross sections. The agreement between theory and 
experiment for the one-jet data is somewhat better than in the TOPAZ case. The low and 
high px data points agree quite well. Only in the intermediate range 3 < px < 5 GeV the 
measured cross section is smaller than the predicted cross section. The two-jet cross section 
agrees very well with the data for all px- We observe that for both experiments the two-jet data 
agree somewhat better with our prediction than the one-jet data. This shows that the NLO 
cross sections correctly account for the experimental data. Also the GRV photon structure 
function which gives a good fit of existing inelastic e7 scattering data produces the correct 
description of the single resolved and double resolved cross sections. Before stronger conclusions 
can be drawn from the comparison one should investigate several effects that could influence our 
predictions. Most of them would change the cross section noticeable only at small px < 3 GeV. 
Such effects are, for example, corrections to the Weizsacker-Williams approximation. Part of 



these corrections have been studied in [12], including intrinsic px effects caused by the non-zero 
momentum of the virtual photons. Other corrections come from the non-vanishing charm mass 
or non-perturbative intrinsic px effects in the photon structure functions. All these corrections, 



as far as they have been estimated || [12] , are below the experimental errors and most of them 
are relevant only at small px- 



1 These are the k factors for the special cone radius R = 1. For smaller or larger cone radii the k factors 
would be completely different 
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4 Conclusions 



Differential inclusive single- and dijet cross sections da/dpr have been calculated in NLO for 
the direct, single resolved and double resolved component as a function of p T . For the double 
resolved two-jet cross section the NLO corrections are estimated with a k factor taken from the 
inclusive one-jet cross section. The sum of these cross sections are compared to the TOPAZ 
and AMY experimental data. We obtained good agreement between measured data and the 
theoretical predictions. For both experiments the two-jet data seem to agree better with the 
theory, than the one-jet data which may be due to the larger experimental errors of the two-jet 
data. The single and double resolved cross sections are obtained with the GRV photon structure 
function which gives also good overall agreement with all measurements on F^- At large pt the 
direct and the single resolved components are the most important ones, so that only at smaller 
Pt the photon structure function can be constrained by the double resolved component. 
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5 Figure Caption 

Fig. 1: Generic diagrams for jet production in 7-7 collisions: (a) direct production in LO, (b) 
direct production in NLO, (c) single resolved production in LO, (d) single resolved pro- 
duction in NLO and (e) double resolved production in LO. 

Fig. 2: da/dp T with R = 1 as a function of px for direct (dashed), single resolved (dashed- 
dotted) and double resolved (dotted) production in NLO. Full curve is the sum of all 
three components compared to TOPAZ data Q: 

(a) Inclusive one-jet cross section with r\ < 0.7, 6 max = 3.2°. 

(b) Inclusive one-jet cross section with rji, r] 2 < 0.7, 6 max = 3.2°. 

The error bars on the data include both statistical and systematic uncertainties. 

Fig. 3: da/dpT with R = 1 as a function of pr for direct (dashed), single resolved (dashed- 
dotted) and double resolved (dotted) production in NLO. Full curve is the sum of all 
three components compared to AMY data 0: 

(a) Inclusive one-jet cross section with 77 < 1.0, O max — 13°. 

(b) Inclusive one-jet cross section with rji, r] 2 < 1.0, 6 max = 13°. 

The error bars on the data include both statistical and systematic uncertainties. 



8 





9 



g 



Fig. 1 a-e 
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